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Abstract — This paper shows experimentally that standard wireless 
networks which measure received signal strength (RSS) can be used 
to reliably detect human breathing and estimate the breathing rate, 
an application we call "BreathTaking". We show that although an in- 
dividual link cannot reliably detect breathing, the collective spectral 
content of a network of devices reliably indicates the presence and 
rate of breathing. We present a maximum likelihood estimator (MLE) 
of breathing rate, amplitude, and phase, which uses the RSS data 
from many links simultaneously. We show experimental results which 
demonstrate that reliable detection and frequency estimation is possible 
with 30 seconds of data, within 0.3 breaths per minute (bpm) RMS error. 
Use of directional antennas is shown to improve robustness to motion 
near the network. 

1 Introduction 

In this paper, we explore the ability to detect and monitor 
breathing using the changes in received signal strength 
(RSS) measured on many links in a deployed wireless 
network. The ability of a wireless network to make 
measurements that can monitor a person's breathing can 
create new opportunities for improving patient moni- 
toring in health care applications. As one example, post- 
surgical patients can die from respiratory depression and 
airway obstruction, which are unfortunately common 
after surgery due the difficulty of correctly dosing seda- 
tives and pain medications administered to a patient [1J. 
Reliable respiration monitoring is critical to detection of 
these conditions J2), [3), QJ. Breathing monitoring also 
has application in diagnosis and treatment for obstruc- 
tive sleep apnea, in which a person experiences periods 
of low breathing rate or long pauses in breathing while 
sleeping |5]. Finally, breathing monitoring may have 
application in detecting sudden infant death syndrome 
(SIDS), which is one of the largest causes of death in 
infants. Parents with a child at high risk for SIDS may 
wish to use a baby breathing monitor to alert them in 
case their child's breathing becomes depressed or stops. 

We use of measurements of RSS between many pairs 
of wireless devices in a deployed network to non- 
invasively detect and monitor a person's breathing, an 
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application we call BreathTaking. While severe fading in 
mobile radio channels is expected, it is counterintuitive 
that small changes in a person's size as a result of 
their breathing could be detected using measurements 
of RSS. However, in this paper, we demonstrate that 
in an otherwise stationary environment, when we use 
the data collected on many links between static wireless 
devices, breathing monitoring is not only possible, but 
remarkably reliable. 
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Fig. 1. Normalized averaged power spectral density 
(PSD) vs. frequency (Hz). Total PSD is defined as the 
argument of ©. The peak of the PSD plot is at 0.253 
Hz (15.18 bpm), compared to the actual breathing rate of 
0.250 Hz (1 5 bpm), shown as a vertical dashed line. 

Our research on BreathTaking is motivated by exper- 
imental observations. We have observed that when a 
person is standing on or near the line-of-sight (LOS) 
of a static radio link, the RSS can be changed simply 
by the person's inhaling and exhaling. In fact, when we 
measure the actual breathing rate and analyze the power 
spectral density of the RSS link data from the network 
as a whole, as shown in Figure [TJ we see a peak very 
close to the actual breathing rate. 

Our BreathTaking does not provide a direct measure of 
breathing. In contrast to End-Tidal CO2 monitoring, for 
example, we do not measure the gasses exhaled from a 
person's nose. We simply make an observation about the 
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presence (or absence) of a strong frequency component 
in the measured RSS in the human breathing range. 
Adults, at rest, breathe at about 14 breaths per minute 
(bpm) [6], while newborns breathe at 37 bpm [7) . To be 
inclusive, we consider a range of 10 to 40 bpm (0.167 
to 0.667 Hz). Few other objects have cyclic motion with 
periodicity in this range, but if there were such an object 
in the deployment region, it might cause the same type 
of observation in the network RSS data. Context remains 
important to interpret results from the proposed system. 

In this paper, we make the following important con- 
tributions in relation to system design, capabilities, and 
limitations of BreathTaking. First, we develop methods 
to accurately estimate the breathing rate and reliably 
detect breathing of a person in the deployment area 
of a wireless network by considering the RSS mea- 
surements on many links simultaneously We approxi- 
mate the breathing signal to be sinusoidal and use the 
maximum likelihood estimation (MLE) to estimate the 
breathing rate. Second, we perform extensive experimen- 
tal evaluation of BreathTaking in an indoor setting. We 
demonstrate monitoring the breathing of an otherwise 
motionless person, in an hospital room with no other 
person present. With thirty seconds of RSS data in a 
twenty-device network, we demonstrate (1) breathing 
rate estimates with RMSE between 0.1-0.4 breaths per 
minute (bpm); (2) a breathing detector without false 
alarms or missed detections during experiments per- 
formed with devices connected to directional antennas. 
We address the performance as a function of the number 
of devices in the network, relative position of the person 
with respect to the sensors, and actual breathing rate. 

Finally, we quantitatively address the following key 
questions that relate to the capability to use RSS mea- 
surements from static wireless networks to monitor 
breathing for the above-described applications: 

1) What is the benefit of using data from multiple 
links simultaneously, as opposed to from one link? 

2) How accurately can breathing rate be estimated? 

3) How long of a measurement duration is required? 

4) What is the effect of the directionality of the anten- 
nas? 

5) How many devices are required for accurate mon- 
itoring? 

6) Is there information in the phase of the breathing 
signal? 

Breathing monitoring in a wireless network has ap- 
plications and implications besides health care. A search 
and rescue team may arrive at a collapsed building and 
throw transceivers into the rubble, hoping to detect the 
breathing of anyone still alive inside. Police or SWAT 
teams may deploy a network around a building to 
determine if people are inside. On the other hand, the 
ability to measure breathing from a wireless network 
has privacy implications. We have shown previously 
that a network deployed around the external walls of 
a building can detect and track a person who is moving 



or changing position |8|, [9]. If this system can also detect 
a person's breathing, it can also detect people who are 
sitting or laying motionless. 

The rest of the paper is organized as follows. In Section 
[2] we present the approach we have used for Breath- 
Taking. Next, we describe the experimental testbed and 
methodology in Section [3] The results of the experiments 
are presented and discussed in Section [4] Finally, related 
work and conclusions are presented in Sections and [6] 

2 Methods 

In this section, we define the measurements, models, and 
goals of the BreathTaking system. 

2.1 Network 

We assume a network in which received power (also 
called RSS) measurements can be made on many links 
between pairs of wireless devices. We assume these 
measurements can be made often, at regular intervals. 
Specifically, assuming a maximum breathing rate of 40 
breaths per minute, BreathTaking requires RSS measure- 
ments to be made at a rate higher than 4/3 Hz, the 
Nyquist sampling rate. 

We denote yi(i) to be the dBm received power on link 
I measured at time i. Note that link I is an ordered pair 
(ti,ri) of the particular transmitter t\ and receiver n for 
link I. We do not generally require full connectivity of the 
network, and instead, assume that connected links are 
numbered from 1 through L, where L is the total number 
of measured links. We wish to maximize L and thus use 
a wireless sensor network with a mesh topology in our 
experiments, although we do not exclude networks with 
other topologies. 

2.2 Signal Model 

In the absence of any motion in the environment of the 
network, we denote 

Vl(i) =Vi + ei(i), (1) 

where yi is the mean RSS for link I and ei(i) is additive 
noise. We assume that the noise on link I is i.i.d. zero- 
mean Gaussian. We also assume that the noise ei(i) is 
independent on different links I. In the presence of a 
breathing person, we assume that the link RSS has an 
additional sinusoidal term, 

yt(i) =yi + Ai cos(2tt fT s i + fa) + ej(i), (2) 

where Ay fy, and / are the amplitude, phase, and 
frequency, respectively, of the periodic component of the 
RSS signal on link I, and T s is the sampling period. We 
assume that the periodic component due to breathing 
would have the same frequency on all links I, and that 
the sampling period is made to be identical on all links, 
thus we do not use a subscript I for frequency /. Phase 
and amplitude are expected to differ between links. 
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2.3 Framework 

We denote the measured signal on link / as a vector, 

y, = [yi(l),...,yi(N)Y 

where N is the total number of samples, and () T indi- 
cates the vector transpose. Because the sampling period 
is T s , the measurement vector corresponds to what we 
call the observation period T, 

T = NT S . 

The observation period T is related to the latency of 
breathing monitoring, since we measure RSS for a dura- 
tion T before obtaining estimates of breathing rate and 
detecting whether or not breathing occurred. 

Our objective from the measurements yi, for / = 
1, ... L, is to detect whether or not a person is breathing 
within the network, and to estimate important parame- 
ters of the model, which we denote 9, 

e = [A T 7 4> T J} T 

where A = [A u . . . , A L ] T and cf> = [</>i, . . . , 4> L ] T . Of 
most interest is the frequency /, as human breathing has 
a characteristic frequency range, from f min to f ma x, as 
discussed in the Introduction. In this paper, we detect 
breathing only within f min — 0.167 Hz to f max = 0.667 
Hz, and in general, we assume that the range is given. 

2.4 DC Removal Filtering 

When estimating the power spectral density of noisy, 
finite-duration y; signals, the mean values yi (the DC 
component) can "hide" the power of lower-amplitude 
sinusoidal components. Yet the DC component does 
not hold information about the presence or absence of 
breathing. Since we assume that breathing is not present 
below a frequency of f m i n , we address this problem sim- 
ply by using a high pass filter which strongly attenuates 
the DC component. We do not require a linear-phase 
filter, but we do need a nearly flat amplitude gain above 
fmin, since we do not want our system to bias towards 
some frequencies because they have been amplified by a 
filter ripple. In all results in this paper, we use a 7th order 
Chebychev high-pass filter with maximum passband 
ripple of 0.1 dB and passband frequency f m in = 0.167 
Hz. An order of 7 was found to be sufficient to have at 
least 40 dB attenuation at frequencies lower than 0.1 Hz. 
For all further discussion, we assume that each link's 
yi{i) signal has been filtered using this high-pass filter, 
and thus yi = 0. 

2.5 Breathing Estimation 

Breathing frequency estimation plays a primary role 
in breathing detection, and thus we discuss frequency 
estimation prior to breathing detection. In this section, 
we present the maximum likelihood estimate (MLE) of 
breathing parameters, including frequency, link ampli- 
tudes, and link phases, given the model presented in 
Section E31 



Maximum likelihood estimation of 8 is an extension of 
the standard sinusoid parameter estimation problem flOl 
p. 193-195] in which there is a single signal composed 
of one sinusoid of unknown phase, amplitude, and 
frequency in additive white Gaussian noise. In our case, 
we additionally have L different link signals, each with 
its own amplitude and phase, and we have a frequency 
limited to the range [f min , f max ] . 

In our case, under the i.i.d. Gaussian noise model 
in the presence of breathing, the likelihood function is 
maximized when the following function J is minimized: 



N 



J W = J2 J2 M*) - A i cos(2ir fT s i + 00] 



(3) 



1 = 1 ! = 1 



One can modify the derivation of flOl p. 193-195] for 
this case, that is, to minimize J (9) in (0, and show that 
a good approximation of the MLE of frequency / is given 
by 

2 

-j2irfTJ 



argmax 

/m«n</</ma* J =1 



N 



(4) 



The approximation is very good whenever the normal- 
ized frequency, fT s , is not very close to or 1/2. In our 
case, we specifically exclude frequencies close to zero, 
and sample at a frequency significantly higher than the 
Nyquist rate. 

Note that if one wishes to estimate breathing fre- 
quency from one link alone, one may use (|D with L = 1 . 

The maximum likelihood link amplitude estimates 
{Ai} and phase estimates {4>{\ are then estimated using 
/, and are given by 



A, 



2 

N 



N 



-j2* fj \i 



(5) 



arctan 



2.6 Breathing Detection 

We consider deciding between two hypotheses: 

Hq : A breathing person is not present (6) 
H-y : A breathing person is present (7) 

For detection, we study two methods: 

1) Single-link breathing detection: Use solely the RSS 
measured on one link in order to detect breathing. 

2) Network-wide breathing detection: Use the RSS mea- 
sured on all L > 1 links in the wireless network to 
detect breathing. 

By comparing the two methods, we quantify the im- 
provement in breathing detection possible when data 
from many links in a network are used. 
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2.6. 1 Single-link breathing detection 

Consider one link's RSS data. Without loss of generality, 
assume the link number is Z = 1. Then the MLE of / 
and Ai are calculated from ||4) and (|5) using L = 1. Our 
simple assumption is that A\ will have higher amplitude 
when a breathing person is present. Thus we detect 
breathing via the hypothesis test, 



Hi 

NA 2 > 



l < Hi 
H 



(8) 



where Ju n k is a user-defined threshold and N is the total 
number of samples. 

2.6.2 Network-wide breathing detection 

From all L measured links, we must decide between Hq 
and Hi . We do not have a statistical model for Ai for the 
case of H\, but we have assumed that the values of Ai are 
higher during breathing. As a first proof of concept, we 
study breathing detection based on a normalized sum of 
the squared amplitudes Af over all links, 



N 



Hi 

< Inet 
1 = 1 H„ 



(9) 



where j ne t is a user-defined threshold and we call S the 
network-wide breathing statistic. Note that S is just a scaled 
version of the maximum of the sum in ((4). Multiplication 
of the average squared link magnitude by N helps 
ensure a constant threshold as a function of N, since 
the average squared link magnitude is approximately 
inversely proportional to N under Hq. 

2.7 Performance Analysis 

We study the experimental performance of each de- 
tector via the probability of false alarm, Pfa, an d the 
probability of missed detection, Pm- The value Pfa is 
the fraction of experiments that do not have breathing 
occurring for which Hi is decided; and Pm is the fraction 
of experiments that do have a person breathing in the 
network for which Hq is decided. Obviously, we'd like 
Pm = Pfa = 0, but there is a tradeoff between the two. 

Clearly, each application will have different require- 
ments for Pm and Pfa- For example, in post-surgical 
breathing monitoring, we will want to have a very low 
Pm, as we do not want to miss the fact that a patient 
has stopped breathing. In contrast, in a search & rescue 
operation, we might be very sensitive to Pfa, because 
saying that there is someone alive in a pile of building 
rubble when there is not would cause us to begin a 
long and fruitless search, when time should be used 
elsewhere to find true casualties. 

We study the performance of breathing rate estimation 
via the RMSE of /, that is, 



RMSE = 



\ 



k=l 



where f(k) and f(k) are the frequency estimate and 
actual breathing frequency, respectively, during experi- 
ment realization k, and there are K total experimental 
realizations. 

3 Experiments 

Our experiments are designed to test BreathTaking for 
use in a medical environment to monitor the breathing 
of a sleeping patient. This section describes the net- 
work hardware and software, environment, experimen- 
tal setup, and actual breathing rate. 

3.1 Network 

We use a network of twenty MEMSIC TelosB wireless 
sensors operating the IEEE 802.15.4 protocol on channel 
26 (a center frequency of 2480 MHz). The sensors run 
TinyOS and SPIN [TT|. a token passing protocol in which 
each node transmits in sequence. When not transmitting, 
nodes are in receive mode, and record the RSS and node 
id for any received packet. Each transmitted packet in- 
cludes the most recent RSS value recorded for each other 
node. A base station sensor, placed in the hallway about 
3 meters from the clinical room door, overhears all of 
the transmitted packet data, which is then time-stamped 
and recorded on a laptop. A packet is transmitted by 
some node approximately once every 12 ms, and thus 
an individual node transmits once every 240 ms. Thus 
each link has its RSS measured at a sampling rate of 4.16 
Hz. We note this is more than sufficient to measure our 
maximum breathing rate of f max = 0.667 Hz. 

To study the effect of antenna directionality, exper- 
iments are conducted with the wireless sensors using 
one of two different antennas: (1) a dipole antenna with 
omni-directional horizontal pattern with gain of 2.25 dBi; 
and (2) a directional patch antenna with a gain of 8.0 dBi. 

3.2 Environment and Setup 

We deploy the network in a clinical room in the Univer- 
sity of Utah School of Medicine. The clinical room is used 
for studies in the Department of Anesthesiology, and is 
designed to appear like a standard hospital patient room, 
with cabinets for blankets and medical supplies, mon- 
itoring equipment, computer monitors, and a hospital 
bed in the center of the room. The hospital bed is a Hill- 
Rom P1900 bed which automatically changes pressure 
in different parts of the bed every 3-4 minutes (and thus 
prevents bed sores in immobile patients). We mention 
this bed movement "feature" because it means that the 
person lying in the bed is not perfectly stationary, even 
when he does not move himself. 

A diagram of the experimental setup is shown in 
Figure |U Sensors are placed on the sides of the beds, 
but not connected in any way to the bed. In this way, 
we ensure that the patient's breathing in no way moves 
the sensors. This is required because we want to show 
that the periodicity in the RSS is caused by the breathing 
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of the person, not the movement of the sensors. Several 
of the sensors are placed at height of 0.9 meters on 
adjustable-height tables positioned at the long sides of 
the hospital bed. Sensors are also placed at a height of 0.8 
meters on two wheeled carts, one metal and one plastic, 
placed at the foot and head of the bed, respectively. 
Finally, four sensors are attached to PVC-pipe stands 
which hold the sensors at a height of 0.9 meters, and 
placed on the sides of the bed. 

3.3 Breathing Rate 

During each experiment in which a person is present in 
the bed, we must have ground truth knowledge of the 
person's breathing rate. In each experiment, the person 
listens to a metronome set to a desired breathing rate, 
and ensures that they breath at the same rate as the 
metronome. The person is also connected to a end-tidal 
CO2 monitor, which involves tubes, two to feed oxygen 
into the person's nostrils, and another two to connect 
the first tubes to a gas sensor which measures CO2 and 
displays it on a screen. The person breathes through their 
nose in all experiments to ensure proper functioning 
of the end-tidal CO2 monitor. The monitor estimates 
frequency from the CO2 sensor's signal, and displays its 
estimated breathing rate. We video record this screen, 
and in post-processing, ensure that the person's breath- 
ing was in fact at the desired rate. 

4 Results 

4.1 Single Link Breathing Monitoring 

In this section, we quantify the observation made in the 
introduction that detection of breathing on any single 
link is unreliable. We compare an experiment run with 
the patch antennas on the nodes, with a person lying 
in the bed with their chest at a height of 1.01 meters 
(Hi) and without a person in the room (Ho). During 
the Hi condition (person lying on the bed) the person 
is using a metronome to breath at a known rate of 15 
bpm. In single-link breathing monitoring, we first use l(4j 
with L = 1 to estimate frequency / which represents the 
breathing rate estimated for an individual link. Then, we 
use © to calculate Ai, which represents the RSS signal 
amplitude at that breathing rate. 

Each link is considered separately, and we estimate Ai 
for each T = 30s observation period for the course of 
the experiments. We use all the links measured during 
the experiment to obtain an ensemble of results that 
characterize the Ai during single-link measurements. 

We display the histograms of Ai in Figure [2] Occur- 
rences are normalized by the total number of realizations 
of Ai in each experiment, and shown on a log scale 
to emphasize the tail behavior. From the histograms, it 
is possible to see that Ai has a heavier tail during Hi 
compared to Hq. However, the maximum Ai recorded 
during Hq is 0.72, and only a very small percentage (63 
out of 9120) of realizations during Hi fall above 0.72. 



That is, a few links, during a few 30-second periods, 
measure higher amplitude sinusoidal components when 
a person is present in the bed than when no one is 
present. If we set a threshold ji mk = 0.72 in (|8) so that 
we have zero false alarms, we would detect breathing 
only on 0.7% of links. 
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Fig. 2. Normalized histogram, on log-scale, of M for 
single link breathing monitoring given H and Hi. Single- 
link amplitudes have a somewhat heavier tail than the no- 
person case. 

Despite the low detection rate, do these links' data 
accurately estimate the breathing frequency? Of the 63 
realizations with A\ > 0.72, 17 provide / estimates 
within 1 bpm of the actual breathing rate (15 bpm). 
However, the other 46 estimates are nearly uniformly 
distributed on [f min , fmax]- 

In sum, one cannot expect to detect breathing or 
estimate breathing rate on any single deployed link. 
Further, even if one has many links in a network, it 
would be unreliable, as a monitoring method, to look 
at single-link amplitude and frequency estimates. 

4.2 BreathTaking Rate Estimation 

Next, we consider at network-wide breathing monitor- 
ing. For the same Hi experiment described in Section l4~Tl 
we study network-wide breathing rate estimates, which 
are calculated using ((4) with L set to the total number 
of links in the network. We calculate rate estimation 
performance for a variety of different periods T. The 
vast majority of breathing rate estimates fall within 5 
bpm of the actual rate - a small fraction do not. We call 
these estimates that are more than 5 bpm from the actual 
rate "invalid" rate estimates, and report the percentage 
of rate estimates that are invalid. We also calculate the 
RMSE and bias of the estimates that are valid, i.e., within 
5 bpm of the true rate. 

The experimental results, shown in Figure [3j show 
that for T > 30s, less than 2% of rate estimates can 
be described as invalid. The RMSE for valid estimates 
is lower than 0.5 for all observation periods T > 25s. 
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Fig. 3. RMSE and bias of valid BreathTaking frequency 
estimates for controlled breathing / patch antenna exper- 
iment (actual breathing rate of 15 bpm), and percent of 
rate estimates that are invalid (right y-axis), vs. observa- 
tion period T. 



1.5 



" 0.5 



-0.5 




0.5 1 
X Coordinate (m) 



— Bed Outline 

— Person 
,5 Sensor 

Coord & No. 

HighAmpl. 

Link 



Fig. 4. Experimental layout showing sensors, bed, and 
person's approximate position. Dashed lines indicate high 
amplitude links. 



For perspective, current medical devices typically re- 
port breathing rate as an integer number of breaths 
per minute, thus rate errors less than 0.5 bpm would 
be insignificant. The bias is small, on the order of 0.1 
bpm. For T > 50s, there are no invalid breathing rate 
estimates. 



4.3 BreathTaking Amplitude Estimation 

Once we obtain the MLE of frequency / using the 
network RSS data as in Section l4~2l we can estimate the 
amplitudes Ai . Note that there are no "actual" values of 
A\) some links will measure high amplitude, and others 
will not. We are particularly interested in the links I 
with particularly high A\. For the same Hi experiment 
described in Section 14.11 consider the Ai for T = 30 s. 
Only 5% of links I have an amplitude Ai > 0.331, links 
we refer to as high amplitude links. 

In Figure 3J we plot the locations of the high amplitude 
links by drawing a dashed line between their transmitter 
and receiver coordinates. We can see that the links 
which cross through the chest area are the ones that are 
particularly affected by breathing. 

Still, only a fraction of the links that cross through 
the chest measure high A{. What other characteristics 
do these high amplitude links have? We find that high 
amplitude links have unusually low average RSS. Over 
all links, the average measured RSS is -39.9. When 
considering only high amplitude links, the average mea- 
sured RSS is -48.6, almost 9 dB lower, a very significant 
difference. This difference cannot be explained by longer 
path length - high amplitude links are only 13% longer, 
on average, than the average path length of all links. 
Since links that in a deep fade experience greater tem- 
poral variations due to changes in the environment (9), it 
makes sense that the amplitude of the breathing-induced 



change would be more noticeable for links with lower 
than average RSS. 

One lesson is that links in a deep fade are more use- 
ful for breathing monitoring. Future work may explore 
changing the center frequency on links, or measuring 
wideband frequency response, with the goal of adap- 
tively using data from links' frequency nulls. Although 
we do not explore this idea in this paper, we would 
expect to be able to improve results by taking advantage 
of "deep fades" wherever in the frequency spectrum they 
occur. 



4.4 BreathTaking Detection Performance 

In the network-wide case, breathing detection is per- 
formed using a normalized sum of the squared am- 
plitudes Ai | 2 over all links I, as given in ((§)• In this 
section, we study the performance of this detector. Using 
the same experiments as presented in Section 14.11 we 
calculate S from (|9]l for each T second period, testing 
detector performance for each T in the range 15 to 60 
seconds, in 5 second increments. Figure |5] shows the 
probability density functions (pdfs) of S for H and Hi 
cases, for T = 15 (top subplot) and T = 30 (bottom 
subplot). 

We find first that S in the Hq case always falls in a 
narrow range, between 0.98 and 1.45. During Hi, the 
S values has a minimum of 1.57 (at T = 15). Further, 
as T increases, S values also increase. For T = 30 and 
T = 60, the minimum S values recorded are 1.63 and 
2.03, respectively. 

We can conclude that for this experiment, because of 
the lack of overlap in value of S in the two cases, one 
can build a reliable detector. For example, we can set 
7„ e t = 1.50 in 10 and perfectly distinguish between no 
person present and a breathing person present. 
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T = 15s Network-wide Breathing Statistic S 
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Fig. 5. Probability density functions for S given H 
(Empty) vs. given Hi (Person) for the (top) T = 15 
and (bottom) T = 30 second observation periods, for 
experiments using patch antennas. 
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Fig. 6. Probability density functions for the S given H 
(Empty) vs. given Hi (Person) for the (top) T = 15 
and (bottom) T = 30 second observation periods, for 
experiments using dipole antennas. 



4.5 Antennas 

Experiments discussed above use patch antennas with 
a 8 dBi antenna gain. In this section, we compare the 
results when using less directional antennas. We run 
more experiments with all nodes connected to dipole 
antennas, which have a omnidirectional pattern in the 
horizontal plane. The first experiment is an Ho experi- 
ment, with no person in the room, and little movement in 
the hallway outside of the room. The second experiment 
is an Hi experiment, with the person lying in the bed 
and breathing at a constant rate of 15 bpm (using a 
metronome). 

The results, shown in Figure [6j are worse than with 
the directional antennas. The values of S given Hq have 
increased regardless of T, to the range from 1.15 to 1.85. 
The values of S given Hi now overlap with those given 
Hq for T = 15s (and T up to 25 s), so regardless of the 
threshold chosen, we cannot have a perfect detector. 

We further believe that movement in the hallway 
outside of the room will have a greater impact when 
using dipole antennas, as compared to when using patch 
antennas. To test this, we have an experimenter stand 
outside of the (closed) door waving his arms above his 
head and moving from side to side. While this motion is 
occurring, we run two more Hq experiments, one with 
dipole antennas, and one with patch antennas. Using this 
data, we recalculate the values of S given Ho. The results 
are shown in Figure [7] 

The change in Hq values for S changes slightly when 
using patch antennas, and changes significantly with 
dipole antennas. With patch antennas, the maximum 
value of S given H has increased to 1.52 (compared 
to 1.45 without motion outside of the door). With dipole 
antennas, the max has increased to 2.08 (compared to 
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Fig. 7. Probability density functions for the S when 
T = 30s given H (Empty) with motion outside of the door 
vs. given Hi (Person), for experiments using (top) patch 
antennas and (bottom) dipole antennas. 



1.80 without motion outside of the door). In this latter 
case, even for medium T (T = 30s is shown in Figure 0, 
there is overlap in the pdfs of S given Hq and Hi, and 
thus it is impossible to select a threshold for zero false 
alarms and zero missed detections. 

In conclusion, it remains possible to monitor breathing 
using dipole antennas. However, it is clearly better, 
in terms of robustness to motion occurring outside of 
the room, to use directional antennas. Our experiments 
using patch antennas are only marginally affected by this 
motion "noise". 
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Actual Breathing Rate (bpm) 


12.0 


15.0 


19.0 


Percent Estimates Invalid 


0% 


1.6% 


0% 


RMSE of Valid Ests. (bpm) 


0.08 


0.42 


0.30 


Bias of Valid Ests. (bpm) 


-0.03 


-0.16 


0.21 



TABLE 1 

Rate estimation performance for three breathing rates 



Person chest height (m) 


0.88 


1.01 


1.13 


Percent Estimates Invalid 


0% 


1.6% 


0% 


RMSE of Valid Ests. (bpm) 


0.17 


0.42 


0.32 


Bias of Valid Ests. (bpm) 


0.008 


-0.16 


0.24 



TABLE 2 

Rate estimation performance for three bed heights 
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Number of Network Nodes 



4.6 Rate Changes 

In this section, we compare frequency estimation perfor- 
mance when the actual breathing rate changes. First, we 
perform three experiments with sensors using patch an- 
tennas in which the person's breathing rate is either 12.0, 
15.0, or 19.0 bpm. Using the metronome set at different 
rates, the person ensures that their breathing follows the 
desired breathing rate. We set T = 30s for all results. The 
results in Table [l] show best performance at the lowest 
breathing rate tested, 12 bpm, but performance does not 
strictly degrade with increased actual breathing rate. 

4.7 Bed Height 

In this section, we analyze three experiments with the 
bed (and thus the person) at different heights. In the 
experiments discussed to this point, the person's chest 
is at a height of 1.01 m above the ground. Here, we raise 
or lower the bed height so that the person's chest is at 
0.88 m, 1.01m, and 1.13m, in three different experiments. 
At the lowest height, the sensors are nearly line-of- 
sight (LOS), that is the line connecting two sensors are 
mostly unobstructed by the person's body. At the highest 
height, the sensors are predominantly at the height of the 
mattress. In all experiments, the actual breathing rate 
is 15 bpm, and we use T = 30s. Table shows the 
breathing rate estimation performance. We find that the 
best performance is in the nearly LOS case, about half the 
RMSE of the 2nd best height. Note that all bed heights 
show acceptable results, with a very small chance of 
invalid estimate, and RMSE below 0.5 bpm. 

4.8 Fewer Sensors 

In this section, we analyze what happens when we use 
a smaller number of sensors in the network. Since the 
number of links is proportional to the square of the 
number of nodes, we expect that having more sensors 
will dramatically improve performance. In fact, the two 
sensor case is a limiting case which we have explored in 
Section |4~TI which showed that one link is insufficient to 
reliably detect and monitor breathing. 



Fig. 8. RMSE and bias of valid BreathTaking frequency 
estimates for controlled breathing / patch antenna experi- 
ment (actual breathing rate of 15 bpm and person height 
of 0.88 m), and percent of rate estimates that are invalid 
(right y-axis), vs. number of sensors in the network T. 
Results are an average over 100 trials using randomly 
selected subsets of the given size. 

We use the same data collected with the 20-node 
network and test what would have happened with a 
smaller network as follows. Let y c {0, . . . , 19} be the 
subset of nodes which we use in a particular trial. We 
run tests with \y\ = 7, 10, 13, 16, 19. For each subset size, 
we run 100 trials and average the results. Subsets are 
randomly selected in each trial from the full set of nodes 
{0,...,19}. 

Figure |8] presents the results for the RMSE, bias, and 
percent valid of the breathing rate estimates. We can see 
that estimates with just seven sensors in the network 
are poor - almost one in four breathing rate estimates is 
invalid (greater than 5 bpm error). When the number of 
sensors is increased to thirteen, the percentage of invalid 
estimates is 1.3%, and the RMSE of valid estimates is 
below 0.3 bpm; and the results improve somewhat as 
the network size increases to 16 and 19. Notably, there 
are zero invalid frequency estimates with 19 nodes. 

4.9 Phase Estimation 

Beyond link amplitudes and network-wide frequency 
estimates, there is also information to be gathered in 
the phase of the sinusoidal signal due to the person's 
breathing, in particular, for links which have a high 
amplitude \Ai\. Consider that if two links measure a 
sinusoid caused by one person's breathing, both should 
be synchronous, that is, rise or fall at the same times. 
However, we do not know whether inhaling will in- 
crease or decrease the RSS on any particular link, so one 
link may reach a maximum while another link reaches 
its minimum. In terms of phase, the {4>i}i might be it 
radians apart from each other. 

To show this effect in the data, we study the <j>i 
estimates for T = 60s observation period. For example, 
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90° 
5 

4 




270° 



Fig. 9. Estimated phases, fa for high-amplitude links I, 
shown on polar plots for five different experiments. Each 
experiment's fa values, to save space, are plotted on a 
different concentric circle labeled by experiment number 
G {1, 2, 3, 4, 5}. Within an experiment, phases are seen to 
be bimodal, with modes separated by 180°. 

consider what we label experiment 1, a patch antenna 
experiment with the person at height 1.01m breathing 
at 15 bpm. For this experiment, we plot fa for links 
with the highest 5% of amplitudes Ai in Figure [9] in the 
circle with radius 1. One may observe from the figure 
that the phases are bimodal with two modes at about 
90 and 270 degrees. We repeat this plotting for four 
other experiments, denoted experiments 2 through 5, on 
different concentric circles on the same polar plot to save 
space. Again, for each experiment, phase estimates for 
the highest 5% amplitude links have two modes, in each 
case, separated by 180°. 

These results are promising in the sense that when 
there are multiple stationary people breathing in the 
network, we might be able to estimate the number of 
people, even when the breathing rates are nearly identi- 
cal, by the number of modes in the distribution of phase. 
Further, we may be able to improve breathing detection 
and monitoring using the fact that phase distribution is 
bimodal. Such methods must be explored in future work. 

5 Related Work 

Breathing monitoring via capnography is standard prac- 
tice for anaesthetized patients in emergency departments 
and in intensive care units [12). The capnometer uses 
an infrared absorption gas analyzer to measure the 
carbon dioxide concentration in exhaled air (end-tidal 
CO2). The breathing rate in this case is determined by 
finding the frequency content of the CO2 concentration 
signal. The exhaled air is measured using tubes in the 
nostrils {e.g., nasal cannula), or from tubes connected to 
a face mask or tracheal tube, which are connected to the 



capnometer. These tubes may become detached; if so, the 
capnometer will detect apnea and alarm. Generally, the 
mask or nasal cannula may be uncomfortable and limit 
the patient's movement. We propose a new non-invasive 
sensor {i.e., sensor not physically attached to the patient) 
for repiratory monitoring, which would allow a patient 
to sleep normally while being monitored. We note that 
capnography directly measures exhalation, while our 
method indirectly measures breathing via the periodic 
changes in the patient's body size. 

Breathing monitoring can also be performed us- 
ing plethysmography (respiratory inductive or thoracic 
impedance plethysmography). These methods measure, 
using electrodes placed on the body, the change in 
inductance or impedance caused by inhalation and ex- 
halation. These electrodes can be contained in a band 
worn around the chest. This is a method used in home 
monitors for infants at risk of SIDS (13| . In comparison, 
the proposed system does not need to be attached to a 
person's body or have wires connected to the person. 

Note that at physician's offices where procedures 
requiring sedation are performed, capnography and 
plethysmography are not typically used due to equip- 
ment costs. Instead, patients are monitored by a pulse 
oximeter, which measures oxygen saturation in blood. If 
a patient stops breathing, oxygen saturation decreases; 
however, the pulse oximeter will detect this desaturation 
only minutes after breathing ceases fl4l . 

Most closely related to the proposed system are other 
proposed non-invasive breathing sensors. Microwave 
Doppler radar systems have been proposed for breathing 
rate estimation [TBI [16). Ultra-wideband (UWB) radar 
has also been proposed for unobtrusive monitoring of 
patient's vital signs ATI ITS) . In fact, UWB radars may 
even be sensitive enough to be able to detect a stationary 
person's heart rate fl9) . In comparison, our proposed 
system uses off-the-shelf wireless devices which are 
significantly lower in cost than radar devices. Rather 
than using a single high-capability radar transceiver, our 
system uses a network of many simple transceivers. 

6 Conclusions 

This paper presents a non-invasive respiration moni- 
toring technique called BreathTaking which uses signal 
strength measurements between many pairs of wireless 
devices to monitor breathing of an otherwise stationary 
person. We present a maximum likelihood estimator to 
estimate breathing parameters, including breathing rate, 
using all of the measured links' RSS data simultane- 
ously. We present detection algorithms based on those 
estimated parameters, and an experimental testbed and 
procedure to validate BreathTaking. 

Using extensive experimental data collected with a 
person lying in a hospital bed, we demonstrate the 
performance of BreathTaking. We find breathing rate 
can be estimated within 0.1 to 0.4 bpm error using 30 
seconds of measurements. We show that the links most 
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affected by breathing are the ones which receive low 
average RSS. Breathing detection is demonstrated to re- 
liably distinguish between the breathing and its absence 
using 15 seconds of RSS data (in patch antenna experi- 
ments) and using 30 seconds of data (in dipole antenna 
experiments), without false alarm or missed detection. 
Antenna directionality is useful to increase robustness 
to external motion. Interestingly, the estimated phases of 
links which are affected by breathing distribution have 
a bimodal distribution with the two modes separated by 
180 degrees. 

If BreathTaking is to be used for medical purposes, 
extensive evaluation on many people, and in many 
settings, must be performed. In addition, this work 
explored RSS-based breathing monitoring using 2.4 
GHz 802.15.4 radios - the system may benefit from 
transceivers with other physical layer protocols and 
center frequencies. Finally, this sensor may be only one 
sensor used in a monitoring system, and its use with 
other sensors and sources of context information should 
be explored. 
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